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Abstract — An optimal signal processing algorithm is derived 
for estimating the time delay and amplitude of each seatterer 
reflection using a frequency-stepped CW system. The channel 
is assumed to be composed of abrupt changes in the reflection 
coefficient profile. The optimization technique is intended to 
maximize the target range resolution achievable from any set 
of frequency -stepped CW radar measurements made in such an 
environment. The algorithm is composed of an iterative two-step 
procedure. First, the amplitudes of the echoes are optimized by 
solving an overdetermined least squares set of equations. Then, 
a nonlinear objective function is scanned in an organized fashion 
to find its global minimum. The result is a set of echo strengths 
and time delay estimates. Although this paper addresses the 
specific problem of resolving the time delay between the first 
two echoes, the derivation is general in the number of echoes. 
Performance of the optimization approach is illustrated using 
measured data obtained from an HP-8510 network analyzer. It is 
demonstrated that the optimization approach offers a significant 
resolution enhancement over the standard processing approach 
that employs an IFFT. Degradation in the performance of the 
algorithm due to suboptimal model order selection and the effects 
of additive white Gaussion noise are addressed. 


I. Introduction 

T HE physical data measured by a monostatic frequency- 
stepped CW (FSCW) radar system is a sequence of 
complex reflection coefficients as seen from the aperture plane 
to which the system is calibrated. The sequence of measured 
reflection coefficients may be interpreted as samples of the 
radar channel frequency response. Given these measurements, 
the objective is to determine the target range. The standard 
approach 1 1 1 is to perform an IFFT (Inverse Fast Fourier 
Transform) that produces an estimate of the radar channel 
time domain impulse response. Peaks in the impulse response 
correspond to reflections and their time delay corresponds to 
the range to the source of reflection. The resolution of this 
approach is limited by the measurement bandwidth. In certain 
applications, such as re-entry plasma density measurement [21, 
[3], the bandwidth must be kept small. The resolution offered 
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by the IFFT approach is often unsatisfactory and an alternative 
approach is required. Even if there exists no distinct bandwidth 
constraint imposed by the radar channel, maximization of the 
achievable resolution for a particular measurement bandwidth 
is an important problem. 

The FSCW measurement system is an outgrowth of 
Frequency-Modulated CW (FMCW) radar |4| that was well 
established in the 1970’s and applied in several measurement 
environments such as locating buried PVC pipe [5], subsurface 
water [6|, underground tanks, and the level of liquid within 
them 1 7 ] . In these FMCW systems, a sinusoid is transmitted 
whose frequency is swept linearly with time. The differential 
frequency between the transmitted and received signals is 
converted to a range estimate. An FSCW system transmits a 
sequence of sinusoids at different frequencies and measures 
the steady-state amplitude and phase shift induced by the radar 
channel at each discrete frequency. One significant benefit of 
performing the measurements at discrete frequencies is the 
ease with which digital signal processing techniques may be 
applied to the data. With the availability of relatively low-cost 
digitally controlled oscillators, FSCW systems have been the 
focus of much attention in recent years. In 1982, Gjessing et 
ai 1 8 1 presented results combining the data produced by the 
analog FMCW and digital FSCW systems in the context of 
adaptive detection and identification of low-flying aircraft in 
a sea clutter baekround using a matched filter approach. In 
1990, Van Hamme [9 1 presented an FSCW technique based 
upon recursively minimizing a nonlinear objective function 
through approximation using a Taylor series expansion and 
demonstrated its performance through range detection of 
transmission line discontinuities. In 1991, Maricevic et aL 
[ 10] demonstrated that the Matrix Pencil algorithm developed 
by Hua et ai. 1 1 1 1 could be applied to FSCW measurements of 
a transmission line having step discontinuities in characteristic 
impedance with significant resolution enhancement over 
standard IFFT processing. In 1992, Carriere and Moses 
[12] applied a total least squares Prony Estimator to FSCW 
measurements of scale models of commercial aircraft and 
demonstrated a significant enhancement in characterizing radar 
signatures over the Inverse Discrete Fourier Transform. 

This paper presents the derivation of an optimal signal 
processing algorithm for obtaining the parameters of a prop- 
agation channel model. The extracted channel parameters are 
the round-trip time delay and amplitude of each reflection. 
The channel is assumed to be composed of abrupt changes in 
the reflection coefficient profile. The optimization technique 
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is intended to maximize the target range resolution achievable 
from any set of frequency-stepped CW radar measurements 
made in such an environment. The derivation is general in 
the number of echoes, and a constant frequency interval 
between the measurement samples is not required. However, 
in order to illustrate the performance of the optimal processing 
algorithm through its geometric interpretation, the specific 
problem of resolving the time delay between the first two 
echoes is addressed. It wilt be shown that the optimization 
process can be interpreted as scanning a nonlinear objective 
function surface in an organized procedure to find its absolute 
minimum. In the two-reflection problem, the objective func- 
tion surface is 3-dimensional, allowing visual inspection. By 
visually examining the structural change in these 3-D surfaces 
as measurement parameters are changed, significant insight 
into the underlying physical process becomes available at a 
glance. 

The resolution limitations of the IFFT approach are demon- 
strated along with the enhancement offered by the optimal 
processing algorithm using physical measurements made with 
an HP-8510 network analyzer. Several sets of measured data 
are analyzed to illustrate the effects of varying bandwidth, 
target range, model order, and noise level. The first data set 
analyzed is a sequence of 801 reflection coefficients measured 
at the input port of the 7- mm Beatty Standard |13J over 
a bandwidth of 18 GHz. The IFFT performs well in this 
case, clearly resolving two known primary reflections with 
the relatively large measurement bandwidth. The optimization 
approach in one delay dimension is then applied to this same 
data set to illustrate that “correct” model order selection is not 
required for accurately measuring the range to multiple sources 
of reflection. The effects of severe additive white gaussian 
noise are examined. It is shown that noise causes a decrease 
in the depth of the objective function minima as well as a 
perturbation in their location. 

In order to illustrate the resolution enhancement achiev- 
able by optimal processing over the IFFT, the measurement 
bandwidth is reduced from 18 to 2 GHz. With 2-GHz band- 
width, the IFFT approach clearly fails to resolve two distinct 
reflections that are known to exist. The optimal processing 
algorithm is then applied (in two delay dimensions) to the 
reduced bandwidth data and the results reveal the potential 
resolution achievable. The remaining data sets are extracted 
from measurements made at the input port of a test apparatus 
composed of a dielectric interferer and metal plate target, also 
driven by an HP-8510 network analyzer. Results for these data 
sets are similar to those presented for the Beatty Standard. 

This paper is organized as follows. In the next section 
the optimal processing algorithm is derived in detail. The 
algorithm is formulated in the frequency domain and an 
iterative procedure is outlined that decomposes the problem 
into a sequence of linear and nonlinear minimization problems. 
Then, experimental data is presented and analyzed using both 
the standard IFFT and optimization approaches, and their 
performance is compared. Finally, the last section presents 
conclusions that may be drawn from the experimental results, 
and future research directions including unresolved problems 
are presented. 


II. Derivation of the Optimal 
Signal Processing Algorithm 

Given a set of frequency response measurements generated 
by an FSCW radar system, the objective is to determine 
the target range. This may be accomplished by assuming an 
underlying physical model for the channel, computing the 
frequency response of the model at the frequencies at which 
the measurements were taken, and finally minimizing the norm 
of the difference between the physical measurements and the 
values produced by the model. The analytical expression lor 
the objective function to be minimized, denoted by ./, is given 
by 

J =|| II (j (u/’o + zAu>)) — + ?Au>)) || (1) 

where the + /Au)) are the measured complex 

reflection coefficient pairs, the H(j(wo + iAw)) are the values 
of a theoretical model, || • || denotes vector norm (14], and / 
is the index number of each measurement which spans the 
integer range from l — n. J is a function of the parameters ol 
the physical model yet to be specified. 

The physical propagation path that will be assumed contains 
a highly reflective target, possibly embedded in an inhomoge- 
neous medium. The inhomogeneity may be either continuously 
tapered, or composed of distinct media layers each having 
a different permittivity, or any combination of the two. The 
simplest realistic model for the system impulse response has 
the form 

h(t) = - / 1 ) + A 2 t)(t - t 2 ) + Atf(t - + ■ • • (2) 

where the A; are the reflection amplitudes and the f t are the 
time delays to the reflections. The corresponding frequency 
response of the system is obtained by taking the Fourier 
Transform of (2), which produces 

H(ju>) = A x <- juU + A 2 er iuU + + ••• (3) 

Samples of this frequency frequency response form the model 
values comprising + ?'Au;)) in (1). 

+ iAa;)) = .4 1 r- jUu+ ' Au ' ,)<1 4- A 2 e-^" +i ^ )h 
+ + * • • (4) 

Therefore, using the assumed model for the system impulse 
response (2), the objective function ./, (1) is a function of the 
reflection amplitudes and delays. 

Using the 2-norm as the metric in (1), the general least 
squares optimization problem becomes finding the set of 
amplitudes and delays in the model that minimize 

n 

J = ^ " \H(j(uj() -(- /Aw)) — H rn (j(uj o + ?’Au;))| 2 (5) 

i=i 

where the summation index i ranges over the n measured 
data pairs. The objective is to find the global minimum of 
the performance metric (5) with respect to the amplitude-delay 
pairs AiJ'i in the model (2), (3). This objective is achieved by 
decomposing the performance metric (5) into a sum of linear 
and nonlinear components. The linear portion of the objective 
is a set of simultaneous linear equations whose solution yields 
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the optimal amplitudes. The metric (5) is quadratic in the 
amplitude parameters A ( . This allows the A, to be calculated 
in closed form by solving the least squares (LS) problem 

Fa ^ J (6) 

where the complex (uxN) matrix F is composed of frequency 
response estimates based on the model (4) and is given by 
(7), shown at the bottom of the page, n is the number of 
frequency measurements and N is the number of A,. pairs 
to be extracted. The complex vector / is composed of the 
physical frequency response measurements and is given by 

/ = [ff A /(wi) H.u(uj 2 ) H m(ujj) ■ ■ ■ H\i ( )] r (8) 

The vector a in (6) contains the set of amplitudes A , to be 
determined 

a = [A! A 2 -L • • • Tv] 7 (9) 

Solution of the least squares problem requires assuming a set 
of delays t It will be shown that by ‘‘scanning “a sequence of 
amplitude optimized objective function values for certain delay 
combinations, it is possible to find the global minimum of (5). 
Finding the global, or absolute minimum of (5) determines the 
optimal set of reflection amplitude-delay pairs. If there exists 
any a priori knowledge about the target range such as upper 
and/or lower bounds, the range of the time delays over which 
the objective function is scanned can be limited and the amount 
of computation necessary to find the global minimum reduced. 
Placing bounds on the time scan of the objective function 
corresponds to time-gating, a procedure often implemented in 
conventional pulsed radar [16]-|21|. 

The solution to (5) will yield a set of complex amplitudes 
because the model (2-4) will not fit the physical data exactly. 
However, in the derivation which follows, the amplitude vector 
a is constrained to be real. The same approach that is used to 
solve the standard real LS problem using the normal equation 
approach [14|, (15] can be extended to solve this problem. 
The solution to (5) is found by minimizing the real scalar that 
results from the following squared inner product. 

y= || (F r +jF l )a-(f,.+jf i ) ||| (10) 

V = ( F,a - f r ) T (F ,o - f r ) + (F,a - f,) r (F,a - /,) 

(ID 

V =« T [FlF r + FjF,]a - a r \Fj.f, + Fj /,} 

-[pF r +fjF,)a + fJf,. + fJf l (12) 


In order to minimize the gradient of (12) is taken with 
respect to the amplitude vector a and set equal to the zero 
vector. 

'f' = 0 = 2 [FlF r + FjF,}a - 2[F 7 ’/,. + Fj f ,\ 

Ua 

(13) 

The result of the minimization of (10) is a real, square set 
of linear equations which can be solved by using any linear 
equation solving technique. 

\Fj. F r + Fj F ,]a = F'i'fr + Fjf, (14) 

However, in certain cases, the symmetric matrix premultiply- 
ing a is extremely ill-conditioned [14], [15]. In order to see 
when such cases arise let 

H = Fj. F v + Fj Fj (15) 

and examine the case of a two reflection problem for simplicity 
(see (16), shown at the bottom of the page). One obvious case 
that will cause H to be ill-conditioned occurs when the values 
of t\ and t 2 are nearly equal. When the values of t\ and t 2 
are exactly equal, H becomes singular. In order to deal with 
this problem, the singular value decomposition technique [14], 
[22] is used to solve (14). When any singular value is found to 
be less than 0.01, its value is set equal to zero before the final 
computation of a is performed. This extracts the minimum 
norm solution to (14), which is precisely what is sought. 

The original LS optimization problem (5) has been de- 
composed into a linear problem whose solution yields the 
amplitudes A,, and a highly nonlinear problem whose solution, 
it will be shown, yields the delays It is important to 
recognize that the solution of (14) for the amplitudes A, 
requires that a set of delays t { be utilized in forming the 
entries of the matrix H. 

The following optimization procedure may be followed for 
extracting the A t . I , pairs: 

1) Choose the number N of A* , l t pairs to be determined. 

2) Assume a sequence of /, values and solve (14) for the 
optimal A, . 

3) Search for the global minimum of (5) using the ampli- 
tudes found from the solution of (14). 

It will be shown in the next section that the geometric 
interpretation of the LS objective function ,7, subject to the 
solution of (14), is a multiple minima surface with grooves 
aligned with the delay axes. The grooves in this surface can 
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be searched sequentially and the global minimum found with 
little computation. The global minimum corresponds to the 
set of optimum reflection amplitudes and delays and hence to 
target range. The next section presents FSCW measurements 
made using an HP-8510 network analyzer that are analyzed 
using both the standard IFFT processing technique as well as 
the proposed optimization approach. 


III. Experimental Results 

This section presents FSCW measurements made using 
an HP-8510 network analyzer in two different experimental 
setups. First, measurements of the reflections in an air-line 
coaxial cable section, the 7 -mm Beatty Standard, are presented. 
Then, measured radar data from a test apparatus is examined. 
A comparison of the performance of extracting the reflection 
time delays using the standard IFFT approach, as used by 
the HP-8510 network analyzer internal computer, and the 
optimization approach is presented. 

Consider the 7-mm Beatty Standard 1 1 1, [10], [13] termi- 
nated in a matched load and its theoretical bounce diagram 
as shown in Fig. 1. The Beatty Standard presents two step- 
discontinuities in the characteristic impedance of the line. 
These discontinuities are produced by an abrupt change in the 
outer diameter of the inner conductor. The bounce diagram 
shows the theoretical amplitude and delay of each reflection 
for measurement calibrated to the input port of the Beatty Stan- 
dard. Measurements of the reflection coefficient s n were taken 
over two different bandwidths using an HP-8510 network 
analyzer. The first set of measurements was taken over the 
band 45 MHz- 18 GHz using 801 evenly spaced frequencies. 
The resulting IFFT is shown in Fig. 2. A rectangular window 
[23] was utilized to maximize the resolution between the 
reflections. The 18-GHz bandwidth is large enough to clearly 
distinguish the reflections. In addition, the amplitudes and 
delays of the first three reflections are very accurate. 

The second set of measurements was taken over the band 
4-6 GHz using 801 evenly spaced frequencies, and the re- 
sulting IFFT is shown in Fig. 3. The 2-GHz bandwidth is so 
small that the two reflections (83 ps, 583 ps) have merged into 
a single pulse. In this case the IFFT cannot resolve the two 
reflections. In addition, the location of the peak in the impulse 
response is biased [24] away from the correct location of either 
of the two reflections. Hence, when two reflections are so close 
in time that the IFFT cannot resolve them, accuracy is lost as 
well. 

Consider the objective function (5) for two assumed delays, 
one scanned from 0-0.8 ns along one delay dimension and the 
other from 0.06-0.1 ns using the Beatty Standard data and 18 
GHz bandwidth. This function may be interpreted as a surface 
above the “delay” plane as shown in Fig. 4. It is important to 
recognize that every value of J(tiJj) is the result of solving 
(14) using each delay pair in the domain shown. The 
surface shown in Fig. 4 has a single global minimum that 
corresponds the optimum delay pair (^ = 83 ps, t 2 = 583 
ps). The grooves in the objective function extend along the 
delay axes. The deeper groove at 0.083 ns is due to the larger 
amplitude reflection caused by the first discontinuity of the 
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Fig. 2. Inverse EFT of 801 measurements of the reflection coefficient at the 
input port of the 7-mm Beatty Standard shown in Fig. 1. The bandwidth used 
in this measurement is 18 GHz. 


Beatty Standard coaxial line (amplitude = —0.3333). The 
groove at 0.583 ns is due to the slightly smaller amplitude 
reflection from the second discontinuity (amplitude = 0.297). 
The surface has multiple local minima, hence a simple descent 
algorithm will fail, in general, to locate the global minimum. 
One approach that will be demonstrated here is to use a grid 
search algorithm on the objective function. Grid search is a 
method of finding the global minimum by finding the least 
value of J along one delay dimension while holding the other 
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Fig. 3. Inverse FFT of 801 measurements of the reflection coefficient at the 
input port of the 7-mm Beatty Standard shown in Fig. 1. The bandwidth used 
in this measurement is 2 GHz. 



Fig. 4. Amplitude optimized objective function constructed from 801 re- 
flection coefficient measurements of the Beatty Standard (Fig. 1) using a 
bandwidth of 18 GHz. 


delay fixed. Then, once the least value of J is found in that 
dimension, the corresponding delay is held fixed, and the 
search is performed along the other dimension. By performing 
multiple iterations in this fashion, the global minimum is found 
once the delays are no longer updated. The global minimum 
for the 18-GHz case (Fig. 4) was found in two iterations to be 
t x = 0.082 ns, t 2 = 0.588 ns, very close to both the theoretical 
delays as well as the estimates produced by the IFF! approach. 

The optimal processing technique is globally optimal in the 
sense that the minimized objective function (5) has many local 
minima and a search procedure is used to find the absolute 
minimum. This global minimum determines the best set of 
reflection amplitudes and delays given the frequency domain 
data and the model that assumes that the radar channel has 
distinctly reflective scatterers. 


d(t t ,t 2 ) 



t 2 (ns) 


Fig. 5. Amplitude optimized objective function constructed from 801 re- 
flection coefficient measurements of the Beatty Standard (Fig. I) using a 
bandwidth of 2 GHz. 


Without a priori knowledge of the number of scatterers, 
there is no way to know, at least initially, the optimal dimen- 
sion in which to search for the global minimum. The more 
interesting case occurs when the bandwidth is reduced to 2 
GHz and the 1FFT is useless for determining the range to 
either of the first two reflections. Applying the optimal signal 
processing algorithm to the frequency data for this case results 
in the objective function illustrated in Fig. 5. There still exist 
grooves in the objective function, each of which is aligned 
with a delay axis. However, there are several grooves whose 
minima lie between the values 0.083 and 0.583 ns. There does 
exist a single global minimum in this domain, and the delay 
values obtained using two iterations of grid search to find it 
are 0.082 and 0.587 ns. This result alone shows the resolution 
enhancement offered by the optimization approach over the 
standard IFFT as used by the HP-8510 network analyzer. 

It is often the case that the number of scatterers is unknown. 
If the Beatty Standard 18-GHz data is scanned in only one 
delay dimension using the optimization approach, the objective 
function shown in Fig. 6 results. Three cases are presented 
to illustrate the effect of additive white Gaussian noise (zero 
mean), in addition to the effect of optimizing in a dimension 
less than the number of scatterers. When no noise is added 
to the data, the two primary reflections are clearly evident by 
the minima in the objective function. The locations in time of 
these minima correspond, with high accuracy, to the two-way 
travel times to the discontinuities in the Beatty Standard. 

By scanning the data in one delay dimension and observing 
distinct multiple minima in the objective function, it is possible 
to detect the presence and location of multiple scatters, as 
well as estimate the “proper" dimension in which to perform a 
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Fig. 6. Effects of noise on the amplitude optimized objective function in 
one delay dimension. 801 measurements of the reflection coefficient at the 
input port of the 7-mm Beatty Standard using a bandwidth of 18 GHz. 


Fig. 7. Effects of noise on the amplitude optimized objective function in 
one delay dimension. 801 measurements of the reflection coefficient at the 
input port of the 7-mm Beatty Standard using a bandwidth of 18 GHz. The 
objective functions have been shifted to a common maximum value. 


subsequent optimization. Fig. 6 also illustrates the degradation 
in the one-dimension objective function as the SNR is reduced. 
The SNR is defined here to be 

n 

£ | #*/(«,•) I 2 

SNRju = 10 log (17) 

y>(w,) 

/==L 

where n(uJi) is the Gaussian noise sample added to the 
measured reflection coefficient (a;,) 1 at frequency u>j. As 
the noise level is increased, the mean value of the objective 
function increases due to the power contributed by the noise. 
This trend is evidenced by the three plots shown in Fig. 6. The 
upward shift in the objective function as the SNR is reduced 
does not degrade the detectability of the scatterers. Flowever, 
as the noise level is increased, the depth of the minima 
decreases, and hence confidence of detection is reduced. At an 
SNR of -10 dB, the two discontinuities are still detectable, 
and the locations in time of the objective function minima still 
provide an accurate estimate of their location. Fig. 7 presents 
a comparison of the same three objective functions of Fig. 
6 with each function shifted to a common maximum value. 
This figure illustrates the relative depths of the minima as a 
function of noise level for a single noise realization. As the 
noise level is increased, the ratio of the minimum value to the 
first side-lobe amplitude diminishes, and hence the confidence 
of detection diminishes as well. When the noise power is 
increased from zero to ten times the signal power, the depth 
of the primary minima is reduced by approximately an order 
of magnitude. The optimized amplitude as a function of delay 
for the three cases of Figs. 6 and 7 is shown in Fig. 8. The 
amplitude of the perturbations in the optimization process are 
more clearly defined in this figure. The perturbation amplitude 
clearly increases with the noise level, and the deviation of the 
estimated amplitudes from the correct values is biased toward 
an overestimate. It remains the subject of future research 
to quantify the statistical characteristics of the optimization 

1 In the ease of an actual noisy measurement, } contains the noise. 



Fig. 8. Effects of noise on the optimal amplitude of a single impulse as a 
function of time delay. 801 measurements of the reflection coefficient at the 
input port of the 7-mm Beatty Standard using a bandwidth of 18 GHz. 

process as a function of noise, as well as bandwidth and target 
range. 

In order to test the applicability of the optimization approach 
to detecting the range to a highly reflective target in the 
presence of an interfering reflection using FSCW radar data, 
a test fixture was constructed at NASA LaRC. This fixture 
is shown in Fig. 9 along with its physical connection to 
the network analyzer via a rectangular-to-circular waveguide 
transition. The network analyzer is calibrated to the input port 
of the waveguide as shown in the diagram. The dielectric 
interferer is composed of a 0.25-in. -thick layer of silicon 
dioxide, or quartz glass U r = 3.78), and a 0.5-in. -thick layer of 
polyurethane foam (c r = 1.4). A 0.5-in. -thick aluminum plate 
is used for the aperture ground plane. The target is formed by 
the presence of another 0.5-in. -thick aluminum plate parallel 
to the aperture plane of the circular horn antenna. Target range 
is set by moving the mobile target. 

In the first set of measurements using the test apparatus, 
the target was set at a distance of 12 in. from the outer 
surface of the quartz glass. The two primary reflections will be 
due to the high dielectric constant quartz glass and the metal 
plate target. 801 measurements of the reflection coefficient 
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Fig. 9. Test apparatus. 


were made from 14-18 GHz, spanning a total bandw idth of 4 
GHz. The IFFT of the measured radar data is shown in Fig. 
10. The first three reflections are clearly resolved. The time 
delay corresponding to the first reflection, due to the quartz 
glass interferer, is difficult to calculate accurately because the 
circular horn antenna as well as the parallel plate waveguide 
formed by the two aluminum plates stimulate the creation of 
several TE and TM modes. However, the prediction of the time 
delay to the interfering reflection from the quartz glass is not 
the primary interest here. The primary concern is target range 
resolution and accuracy. Detecting accurately the time delay 
from the reflection due to the quartz glass to the metal plate 
target is the objective. The IFFT shown in Fig. 10 indicates 
there is a delay of 1 .97 ns from the quartz glass to the metal 
plate. Using the free space velocity of propagation, the two- 
way distance of 24 in. in air requires 2.032 ns. Since each 
reflection from the quartz glass is assumed to be a single 
impulse, two unresolved impulses resulting from the front 
and back surfaces, an additional delay of 41 ps is included 
to account for the 0.25 in. of travel through the glass. Thus, 
the target range is defined to be the distance from the center 
of the glass to the metal plate. The total theoretical two-way 
travel time is 2.073 ns. A time delay estimate of 1.97 ns 
is produced by the IFFT. Thus, the error in the time delay 
estimate produced by the IFFT using 4 GHz bandwidth is 103 
ps and corresponds to a range estimate error of approximately 
0.608 in. (1.545 cm). 

The objective function surface for 801 measurements over 
a 4-GHz bandwidth is shown in Fig. 1 1. This surface clearly 
shows grooves along the delay axes. The domain of this plot 
was chosen to display the symmetry of the objective surface. 
As shown, there are two global minima one with G = 3.47 
ns, t 2 ~ 5.43, ns and the other with G = 5.43 ns, f 2 = 3.47 
ns, which are exactly the same solution to the problem. The 
deeper groove at 3.47 ns is due to the reflection from the quartz 


lh(t)l 



Fig. 10. Normalized Inverse FFT of 801 measurements of the reflection 
coefficient at the input port of the test apparatus shown in Fig. 9. The 
bandwidth used in this measurement is 4 GHz, and the target range is 12 in.. 

glass and the other groove is due to the weaker target reflection 
at 5.43 ns. Between the two major grooves lie several small 
amplitude grooves, also aligned with the delay axes. The top of 
the surface is relatively smooth and is characteristic of a “well- 
behaved” objective function. If the number of measurements 
is reduced from 801 to 21 while holding the bandwidth and 
target range constant at 4 GHz and 12 in., respectively, the 
objective function surface that results is shown in Fig. 12. 
Although the top of the surface is not as smooth as using 
801 measurements, the small amplitude perturbations do not 
increase the number of grid search iterations required to find 
either of the global minima shown. The two global minima 
are deep and narrow' enough to clearly indicate the presence 
of two significant reflections. This result shows the relative 
immunity of the optimization approach to reducing the size 
of the data set. 

In the next set of measurements, the target was initially set 
at a distance of 12 in. from the outer surface of the quartz 
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Fig. 11. Amplitude optimized objective function constructed from 801 
reflection coefficient measurements of the test apparatus (Fig. 9) using a 
bandwidth of 4 GHz. Target is at a distance of 12 in. from the outer surface 
of the quartz glass. 

glass. One-hundred and one measurements of the reflection 
coefficient were made from 14-14.5 GHz, spanning a total 
bandwidth of only 500 MHz. The IFFT of the measured radar 
data is shown in Fig. 13. The bandwidth was selected such 
that the IFFT has a resolution marginally able to detect the 
presence and correct location of the first two reflections. When 
the target distance is reduced to 2.25 in. from the outer surface 
of the quartz glass, the resulting IFFT is shown in Fig. 14. 
The narrow bandwidth spanned by the measurement sequence 
is so small that several reflections have merged into a single 
lobe of the IFFT. The 2.25 in. distance from the outer surface 
of the quartz glass to the metal plate target corresponds to 
0.381 ns of two-way travel time. Including 41 ps of delay 
within the quartz glass results in a total theoretical two-way 
travel time of 0.422 ns. Since the metal plate target is so 
close to the interferer, there are many reflections between the 
quartz glass and the aluminum plate. These reflections are 
separated by 0.422 ns in time. The main lobe in the IFFT spans 
approximately 2 ns, enough time to “swallow” nearly 3 of 
these multi-bounce reflections. The time delay corresponding 
to the second major peak in the IFFT of Fig. 14 is 4.8 ns and 
has no value in assessing the range to the target. This is the 
best performance the IFFT has to offer for determining the 
target range under the bandwidth constraint of 500 MHz and 
101 frequency measurements. 

The objective function surface resulting from the data set 
whose parameters are 4 GHz, 101 measurements, and target 
distance 12 in. from the quartz glass is shown in Fig. 15, The 
domain scanned has been chosen once again to display the 
symmetry of the objective function surface. The two global 
minima each provide reflection delay estimates of 3.47 ns and 
5.50 ns, resulting in a target delay estimation error of only 43 



t 2 (n«) 

Fig. 12. Amplitude optimized objective function constructed from 21 reflec- 
tion coefficient measurements of the lest fixture (Fig. 9) using a bandwidth 
of 4 GHz. Target is at a distance of 12 in. from the outer surface of the 
quartz glass. 
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Fig. 13. Normalized Inverse FFT of 101 measurements of the reflection 
coefficient at the input port of the test apparatus shown in Fig. 9. The 
bandwidth used in this measurement is 500 MHz. The Target is at a distance 
of 12 in. from the outer surface of the quartz glass. 

ps. If the edge of the surface defined by t t ~ 8 ns is examined 
closely, three distinct troughs in the surface can be observed. 
The first trough is located at a delay value of 3.47 ns and 
corresponds to the strongest reflection that is due to the quartz 
glass. The second trough is located at a delay value of 5.50 ns 
and is due to the metal plate target. The third trough is located 
at a delay value of 7.51 ns and is due to the first multi-bounce 
reflection between the metal plate and quartz glass. Hence, 
although the surface presented is the result of minimizing a 
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Fig. 14. Normalized Inverse FFT of 101 measurements of the reflection 
coefficient at the input port of the test apparatus shown in Fig. 9. The 
bandwidth used in this measurement is 500 MH/, and the target range is 
2.25 in. 

two-delay objective, the surface contains information enabling 
the location of higher-order reflections. If the target is moved 
to a distance of 2.25 in. from the quartz glass while holding the 
bandwidth and number of measurements constant at 500 MHz 
and 101, respectively, the objective function surface shown 
in Fig. 16 results. This objective function surface may appear 
to lack structure. However, the surface has grooves running 
along the delay axes that can be scanned successfully using 
the grid search algorithm. Evidence of these grooves may be 
seen by carefully examining the right half of the surface. The 
ridges in the surface indicate the presence of the grooves. The 
surface also possesses a macroscopic bowl shape, albeit very 
“noisy,” in the domain shown. There are two global minima, 
either of which produces the delay estimates /] = 3.53 ns, 
t >2 — 3.95 ns. These delay values were found by polishing the 
result of a first-pass grid search with an additional grid search. 
In the first pass, the domain scanned was 1-6 ns in each 
delay dimension using 50 ps increments. After 16 iterations 
of grid search, the resulting time delays are t\ = 3.55 ns, 
t ‘2 — 3.95 ns. Then, a second pass was performed over the 
domain 3. 5-3. 6 ns in the first delay dimension and 3. 9-4.0 ns 
in the second delay dimension using 1 ps increments. After 
12 iterations of grid search in the second pass, the resulting 
time delays are t\ = 3.53 ns, t<i = 3.95 ns. The impulse 
response produced by the optimization technique is shown in 
Fig. 17. The most important measurement is the difference in 
time delay between the two reflections. This delay difference 
is found to be t ,2 - t\ — 0.42 ns. The theoretical value was 
calculated to be 0.422 ns. Hence, a delay error of only 2 ps is 
produced by the optimization approach in this case. This 2-ps 
error in the delay estimate corresponds to approximately 0.3 
mm. This result clearly shows the resolution capability of the 
optimization approach. The IFFT for this measurement case 
(Fig. 14) was totally useless for determining the range to any 
reflection while the optimization approach has distinctly and 
accurately resolved both the reflection from the metal plate 
target as well as the reflection from the quartz glass. 

The results for the IFFT and optimization approaches for 
both the Beatty Standard and test apparatus experiments are 
presented in Table I. 



Fig. 15. Amplitude optimized objective function constructed from 101 
reflection cocflicient measurements of the test apparatus (Fig. 9) using a 
bandwidth of 500 MHz. Target is at a distance of 12 in. from the outer 
surface of the quartz glass. 

IV. Conclusion 

This paper has presented the derivation and demonstration 
of an algorithm for optimal processing of radar data produced 
by any FSCW system. The derivation is based on minimiz- 
ing the two-norm of the difference between the sequence 
of measured reflection coefficients and those produced by 
a model that assumes the echoes from the radar channel 
are impulses in the time domain. The algorithm does not 
require a constant interval between the measurement samples. 
A performance comparison between the IFFT and optimal 
processing techniques was presented in the context of two 
physical measurement sets. In the first set, the Beatty standard 
was terminated in a matched load and driven by an HP-8510 
network analyzer. The measurement bandwidth was reduced 
and it was shown that while the IFFT produced meaningless 
reflection range estimates the optimal processing approach was 
able to clearly resolve the reflections with high accuracy. In the 
second set of measurements, a test apparatus also driven by an 
HP-8510 network analyzer was used to illustrate the effects of 
reducing the bandwidth, number of measurements, and target 
range. It was shown that when the number of measurements 
was reduced to 21, the objective function surface developed 
a high-frequency, low-amplitude perturbation. The surface 
global minima, however, were just as sharp and deep as using 
801 measurements. This displayed the relative immunity of 
the optimization approach to reducing the size of the data 
set. When the target range was reduced from 12.125-2.375 
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TABLE I 

Performance Comparison of the IFFT and Optimation Approaches for 
the Beatty Standard and Test Apparatus Experiments. The Target 
Range in the Test Apparatus Experiments is Defined to bf. the Distance 
from the Center of the Quartz Glass to the Metal Plate Target 



Beatty Standard 




Theoretical 

IFFT 

Optimization 

BW = 18 GHz, n=801 measurements 
Time delay to the first reflection: 

Time delay to the second reflection: 

83 ps 
583 pS 

82 ps 
585 ps 

82 ps 
588 ps 

BW = 2 GHz, rv=801 measurements 
Time delay to the first reflection: 

Time delay to the second reflection: 

83 ps 
583 ps 

wa' 

N/A* 

82 ps 
587 ps 


Test Apparatus 



Target Range = 12.125" 

BW = A GHz, n = 801 measurements 

Theoretical 

JEfl 

Optimization 

Time delay from the quartz glass center 
to the metal plate target: 

2 073 ns 

1,97 ns 

1.96 ns 

Target Range = 12.125" 

BW = 500 MHz, n = 101 measurements 



Time delay from the quartz glass center 
to the metal plate target 

2.073 ns 

1.45 ns 

2.03 ns 

Target Range = 2.375" 

BW = 500 MHz, n = 101 measurements 



Time delay from the quartz glass center 
to the metal plate target: 

0.422 ns 

N/A* 

0 420 ns 


*It is not clear from the IFFT results that there are two reflections. Even if 
two reflections are assumed, the time delays to the peaks in the IFFT results 
are erroneous to the point of being meaningless. 


in. while holding the bandwidth constant (500 MHz), the 
range resolution and accuracy deterioration produced by IFFT 
processing was clearly evident. In this case, the IFFT could 
not distinguish the presence of two reflections nor could 
it determine accurately the range to either the interfering 
reflection or the metal plate target. However, the optimal 
processing approach could not only resolve the two reflections, 
but estimate both of their ranges with high accuracy. This 
result demonstrated the resolution enhancement offered by the 
optima] processing approach over IFFT processing. 

It remains to develop the optimization approach in higher 
dimensions. The simple grid search algorithm may be inef- 
ficient for extracting global minima in higher dimensions, 
and a more sophisticated search method may be required 
[25]. One distinguishing characteristic of the optimization 
approach is its geometrical interpretation which is locating 
the absolute minimum of a surface lying above a delay 
plane. By varying measurement parameters and viewing the 
objective function surface, the effect on the solution can be 
understood intuitively through visual inspection. A potential 
advantage of the optimization approach is its ability to in- 
corporate constraints such as time-gating, a technique often 
implemented in conventional pulsed radar [ 1 6]— [2 1 ]. Time- 
gating is accomplished by searching only over certain regions 
in the delay plane. 

Future research will be directed toward quantifying the 
amplitude and frequency of the perturbations that occur in the 
objective function surface as the bandwidth, or the number of 
measurements, is reduced. The effects of added noise on the 
estimated amplitude and delay parameters will be quantified 



Fig. 16. Amplitude optimized objective function constructed from 101 
reflection coefficient measurements of the test apparatus (Fig. 9) using a 
bandwidth of 500 MHz. Target is at a distance of 2.25 in. from the outer 
surface of the quart/ glass. 
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Fig. 17. Normalized Impulse response of the test apparatus of Fig. 6 
obtained using grid search (16 iterations, 14 iterations polished) on the 
amplitude-optimized objective function resulting from 101 measurements over 
500-MHz bandwidth. The target is at a distance of 2.25" from the outer 
surface of the quartz glass. 

statistically, along with the confidence of target detection. 
The benefit of irregular sampling in the frequency domain 
will be investigated. The number of delay dimensions will be 
increased and other search procedures investigated for efficient 
global minimum scanning. 
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